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ABSTRACT 



Context. The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains contro- 
versial. White dwarfs have been proposed to explain these results and, hence, to contribute significantly to the mass 
budget of our Galaxy. However, several constraints on the role played by regular carbon-oxygen white dwarfs exist. 
Aims. Massive white dwarfs are thought to be made of a mixture of oxygen and neon. Correspondingly, their cool- 
ing rate is larger than those of typical carbon-oxygen white dwarfs and they fade to invisibility in short timescales. 
Consequently, they constitute a good candidate for explaining the microlensing results. 

Methods. Here, we examine in detail this hypothesis by using the most recent and up-to-date cooling tracks for massive 
white dwarfs and a Monte Carlo simulator which takes into account the most relevant Galactic inputs. 
Results. We find that oxygen-neon white dwarfs cannot account for a substantial fraction of the microlensing depth 
towards the LMC, independently of the adopted initial mass function, although some microlensing events could be due 
to oxygen-neon white dwarfs. 

Conclusions. The white dwarf population contributes at most a 5% to the mass of the Galactic halo. 



Key words, stars: white dwarfs — stars: luminosity function, mass function 
matter — Galaxy: structure — Galaxy: halo 
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1. Introduction 

Several cosmological observations show compelling evi- 
dence that baryons represent only a small fraction of the 
total matter in our Universe and that non-baryonic dark 
matter dominates over baryons. To be specific, in the stan- 
dard cosmological model fl\ — 0.72 and — 0.27, 
whereas J7b — 0.044. Moreover, most of the baryons are 
non-luminous, since fi* ~ 0.005. For the case of our own 
Galaxy it has been found that the virial mass out to 100 
kpc is M « 10 12 Mq while the baryonic mass in the form of 
stars is M± w 7 x 10 10 Mq, which means that for the Milky 
Way, the baryon fraction is at most 8% (Klypin et al. 2007). 
This problem is known as the missing bayon problem — see 
the excellent review of Silk (2007) for a complete, interest- 
ing and recent discussion of this issue — and it is critical 
in our understanding of how the Galaxy (an by extension 
other galaxies) were formed and will ultimately evolve. In 
order to solve this problem, three alternatives can be envis- 
aged: either these baryons are in the outer regions of our 
Galaxy, or, perhaps, they never were present in the pro- 
togalaxy or, finally, they may have been ejected from the 
Milky Way. The most promising explanation and the cur- 
rently favored one is the first of these options. 
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The most likely candidates for building up the bary- 
onic dark matter density are massive baryonic halo ob- 
jects, or MACHOs. It has been suggested that MACHOs 
could be planets [M ~ 10 Mq), brown dwarfs (with 
masses ranging from ~ 0.01 to ~ 0.1 Mq), primordial black 
holes (M k, 1O" 16 M ), molecular clumps (M ~ IMq) 
and old white dwarfs (M ~ 0.6 Mq). White dwarfs are 
specially interesting candidates not only because their in- 
trinsic faintness, but also because in addition to the mass 
of the white dwarf itself their progenitors have to return 
to interstellar medium a sizeable fraction of their original 
mass (~ 2 Mq on average) once the white dwarf is formed. 
Additionally, much expectation has been generated since 
the pioneering proposal of Paczyhski (1986) that MACHOs 
could be found through gravitational microlensing. Since 
then, several groups such as the MACHO (Alcock et al. 
1997, 2000), EROS (Lasserre et al. 2001, Goldman et al. 
2002, Tisserand et al. 2006), OGLE (Udalski et al. 1994), 
MOA (Muraki et al. 1999) and SuperMACHO (Becker et 
al. 2005) teams have monitored millions of stars during sev- 
eral years in both the Large Magellanic Cloud (LMC) and 
the Small Magellanic Cloud (SMC) to search for microlens- 
ing events. Among these searches, it is worth mentioning 
that the MACHO collaboration has succeeded in revealing 
13-17 microlensing events during their 5.7 yr analysis of 
11.9 million stars in the LMC (Alcock et al. 2000). In their 
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analysis they derived an optical depth towards the LMC of 
r = 1.2^0 3 x 10~ 7 for events with durations in the range 
2 < t < 400 days. This value is smaller than that expected 
for a full MACHO halo. In fact, it corresponds to a halo 
fraction 0.08 < / < 0.50 at the 95% confidence level with 
a MACHO mass in the range 0.15 Mq < M < 0.50 Mq, 
depending on the halo model. Despite the fact that only 
a fraction of the dark matter could be in the form of 
MACHOs, there is still a large controversy about the nature 
of the reported microlensing events and to which extent 
they contribute to the mass budget of the dark halo of the 
Galaxy. In fact, a large variety of possible explanations have 
been proposed to explain these microlensing events. For in- 
stance, white dwarfs, brown dwarfs and black holes appear 
as natural candidates, whereas self-lensing by stars of the 
LMC (Sahu 1994, Gyuk, Dalai & Griest 2000) has been pro- 
posed as well. Also, other explanations — like tidal debris or 
a dwarf galaxy toward the LMC (Zhao 1998), a galactic ex- 
tended shroud population of white dwarfs (Gates & Gyuk 
2001), blending effects (Belokurov, Evans & Le Du 2003, 
2004), non-conventional initial mass functions (Adams & 
Laughlin 1996; Chabrier et al. 1996), spatially varying mass 
functions (Kerins & Evans 1998, Rahvar 2005), and other 
explanations (Holopainen et al. 2006) — have been also 
thoroughly discussed during the last years. However, all of 
these proposals have been received with some criticism be- 
cause none of them fully explains the observed microlensing 
results. 

There are as well other observations that are important 
pieces of evidence in this puzzle, such as the results of the 
EROS collaboration or the search for very faint objects in 
the Hubble Deep Field. We briefly summarize them. The 
EROS team has recently presented an analysis of a subsam- 
ple of bright stars from the LMC, minimizing the source 
confunsion and blending effects (Tisserand et al. 2006). 
Their results imply that the optical depth towards the LMC 
is t < 0.36 x 10~ 7 at the 95% confidence level, correspond- 
ing to a fraction of halo mass of less than 7%. This result is 4 
times smaller than that obtained by the MACHO team and, 
consequently, sets a strong upper bound to the contribution 
of MACHOs to the mass budget of the Galactic dark matter 
halo. Nevertheless, the nature of the observed microlensing 
events still remains to be clarified. Also, the Hubble Deep 
Field-South has provided another opportunity to test the 
contribution of white dwarfs to the Galactic dark matter. 
In particular, Kilic et al. (2005) have recently found three 
white dwarf candidates among several faint blue objects 
which exhibit significant proper motion and, thus, are as- 
sumed to belong to the thick-disk or the halo populations. 
If in the end these white dwarfs are spectroscopically con- 
firmed it would imply that white dwarfs can account for 
about 10% of the Galactic dark matter, which would be 
consistent with the results of the EROS team, and with 
previous estimates (Chabrier 2004). 

In a previous paper (Garcfa-Berro et al. 2004) we exten- 
sively analyzed the role played by the carbon-oxygen (CO) 
white dwarf population in several different observational re- 
sults, namely, the reported microlensing events towards the 
Large Magellanic Cloud (Alcock et al. 2000), the results of 
the Hubble Deep Field (Ibata et al. 1999) and the results of 
the EROS experiment (Goldman et al. 2002). We performed 
a thorough study for a wide range of Galactic inputs, in- 
cluding different initial mass functions and halo ages, and 



several density profiles corresponding to different halo mod- 
els. Our main result was that a sizeable fraction of the halo 
dark matter cannot be locked in the form of old hydrogen- 
rich white dwarfs with CO cores. Specifically, we found that 
this fraction should be of the order of 4%, in agreement with 
the standard models of the Galactic halo. However in our 
analysis we disregarded the contribution of massive white 
dwarfs, that is, stars more massive than ~ 1.1 Mq. The 
core of these white dwarfs consists of a mixture of oxygen 
and neon. Since oxygen-neon (ONe) white dwarfs cool con- 
siderably faster than the bulk of CO white dwarfs (Althaus 
et al. 2007) it is reasonable to expect that perhaps some 
of the microlensing events could be due to these elusive 
massive white dwarfs. It is also worth mentioning at this 
point that the MACHO collaboration in their first season 
reported a microlensing event with a duration of 110 days 
towards the galactic bulge (Alcock et al. 1995). For this par- 
ticular event a parallax could be obtained from the shape 
of the light curve, from which a mass of 1.3^ ' 1 Mq was de- 
rived, indicating that the gravitational lens could possibly 
be a massive ONe white dwarf or a neutron star. Moreover, 
studies about the distribution of masses of the white dwarf 
population (Finlcy et al. 1997; Liebcrt et al. 2005) show the 
existence of a narrow sharp peak near 0.6 Mq, with a tail 
extending towards larger masses, with several white dwarfs 
with spectroscopically determined masses within the inter- 
val comprised between 1.0 and 1.2 Mq. 

In this paper we analyze if ONe white dwarfs could be 
responsible for a sizeable fraction of the reported microlens- 
ing events towards the LMC. The paper is organized as 
follows. In Sect. 2 we briefly describe the main ingredients 
of our Monte Carlo code and other basic assumptions and 
procedures necessary to evaluate the microlensing optical 
depth towards the LMC. Section 3 is devoted to describe 
our main results, including the contribution of ONe white 
dwarfs to the halo white dwarf luminosity function and to 
the microlensing optical depth towards the LMC, and we 
compare our results to those of the MACHO and EROS 
teams. In this section we also check if ONe white dwarfs 
could be detected in the Hubble Deep Field South and we 
discuss the contribution of ONe white dwarfs to the bary- 
onic content of the Galaxy. Finally, in Sect. 4 our major 
findings are summarized and we draw our conclusions. 

2. The model 

2.1. The Monte Carlo simulation 

An extensive description of our Monte Carlo simulator 
has been already presented in Garcfa-Berro et al. (2004). 
Consequently, here we will only briefly summarize the main 
ingredients of our model. We have used a random number 
generator algorithm (James 1990) which provides a uniform 
probability density within the interval (0, 1) and ensures a 
repetition period of il 10 18 , which is enough for our pur- 
poses. Each one of the Monte Carlo simulations discussed 
in Sect. 3 below consists of an ensemble of 40 indepen- 
dent realizations of the synthetic white dwarf population, 
for which the average of any observational quantity along 
with its corresponding standard deviation were computed. 
Here the standard deviation means the ensemble mean of 
the sample dispersions for a typical sample. 

We have considered an otherwise typical spherically 
symmetric halo. The density profile of this model is the 
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isothermal sphere of radius 5 kpc, also called the "S- 
model", which has been extensively used by the MACHO 
collaboration (Alcock et al. 2000; Griest 1991). Despite the 
existence of other density profiles, such as the exponen- 
tial power-law model, the Navarro, French & White (1997) 
density profile and others, in our previous study (Garcia- 
Berro et al. 2004) we showed that the differences between 
them are not significant for the case under study and, conse- 
quently, we adopt the most simple description. The position 
of each synthetic star is randomly chosen according to this 
density profile. 

We have considered two different initial mass functions, 
the rather standard initial mass function of Scalo (1998) 
and the biased log-normal initial mass function proposed 
by Adams & Laughlin (1996), which is very similar to the 
non-conventional initial mass function of Chabrier et al. 
(1996). This biased initial mass function has been included 
just for the sake of completeness, since it does not seem 
to be compatible with the observed properties of the halo 
white dwarf population (Isern et al. 1998; Garcfa-Berro et 
al. 2004), with the contribution of thermonuclear super- 
novae to the metallicity of the Galactic halo (Canal et al. 
1997), and with the observations of galactic halos in deep 
galaxy surveys (Chariot & Silk 1995). The main sequence 
mass is obtained by drawing a pseudo-random number ac- 
cording to the adopted IMF. Once the mass of the pro- 
genitor of the white dwarf is known we randomly choose 
the time at which each star was born. We assume that the 
halo was formed 14 Gyr ago in an intense burst of star for- 
mation of duration ~ 1 Gyr. Given the age of the halo, 
the time at which each main-sequence progenitor was born 
and the main sequence lifetime as a function of the mass 
in the main sequence (Iben & Laughlin 1989) we know 
which stars have had time enough to enter in the white 
dwarf cooling track, and given a set of theoretical cooling 
sequences and the initial to final mass relationship (Iben & 
Laughlin 1989), we know their luminosities, effective tem- 
peratures and colors. The cooling sequences adopted here 
depend on the mass of the white dwarf. White dwarfs with 
masses smaller than Mwd = 1-1 Mq are expected to have a 
CO core and, consequently, for them we adopt the cooling 
tracks of Salaris et al. (2000). White dwarfs with masses 
larger than Mwd = 1.1 Mq most probably have ONe cores 
and for these white dwarfs we adopt the most recent cool- 
ing sequences of Althaus et al. (2007). Both sets of cooling 
sequences incorporate the most accurate physical inputs 
for the stellar interior (including neutrinos, crystallization, 
phase separation and Debyc cooling) and reproduce the 
blue turn at low luminosities (Hansen 1998). Also, the en- 
semble of cooling sequences used here encompass the full 
range of interest of white dwarf masses, so a complete cov- 
erage of the effects of the mass spectrum of the white dwarf 
population was taken into account. 

The kinematical properties of the halo white dwarf pop- 
ulation have been modeled according to a gaussian law 
(Binney & Tremaine 1987): 



where ay and a t — the radial and the tangential veloc- 
ity dispersion, respectively — are related by the following 
expression: 
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which reproduces the flat rotation curve of our Galaxy 
at large distances. We have adopted a circular velocity 
V c = 220 km/s. Finally, and in order to obtain the he- 
liocentric velocities we have take into account the peculiar 
velocity of the sun {Uq,Vq,Wq) = (10.0,15.0,8.0) km/s 
(Dchnen & Binney 1998). Since white dwarfs usually do 
not have determinations of the radial component of the ve- 
locity, the radial velocity is eliminated when a comparison 
with the observational data is needed. Moreover, we only 
consider stars with velocities larger than 250 km/s because 
white dwarfs with velocities smaller than this would not 
be considered as halo members. Additionally, we also dis- 
card stars with velocities larger than 750 km/s, because 
they would have velocities exceeding 1.5 times the escape 
velocity. 

2.2. Modeling the microlensing events towards the LMC 

In order to ascertain the contribution of halo white dwarfs 
to the microlensing events towards the LMC we have pro- 
ceeded in three steps. First of all we have built a model of 
the LMC following closely the procedures detailed in Gyuk 
et al. (2000) and Kallivayalil et al. (2006). This model takes 
into account, among other parameters, the scale length and 
scale height of the LMC, its inclination and its kinematical 
properties. This model provides us with a synthetic popula- 
tion of stars representative of the monitored point sources. 
In a second step we search for those halo white dwarfs that 
could be responsible of a microlensing event. This implies 
that the white dwarf should be fainter than a magnitude 
limit, otherwise it would not be considered as a genuine mi- 
crolensing event. Typically we have taken m™' = 17.5 mag , 
which is the value adopted by Alcock et al. (2000). This 
value has been confirmed to be a reasonable estimate by 
the detailed theoretical simulations of Garcfa-Berro et al. 
(2004). Finally, we check if the angular distance between 
the white dwarf and the monitored star is smaller than the 
Einstein radius 9e = Re/Dol, where Dql is the distance 
between the observer and the lens and Re is the Einstein 
radius which is given by the expression 
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where Dos is the observer-source distance and x = 
D 'ol /Dos- If this is the case then we have a microlensing 
event and we compute the corresponding probability. This 
probability is integrated over the total monitoring period 
of observation and filtered by the detection efficiency func- 
tion, which allows us ot obtain the optical depth (Alcock 
et al. 2000): 
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where E is the total exposure (in star-years), ti is the 
. . Einstein ring diameter crossing time, and e(ti) is the de- 
\ ' tection efficiency The detection efficiency and E depend 



4 



J. Camacho et al.: ONe white dwarfs and the Galactic halo 



-2.0 



co 
I 

O 

Oh 



SO 

o 



-4.0 - 



-6.0 - 



o 

J3 



CO 
I 

O 

a 



o 



-6.0 - 



I i i | i i i | i i i | i i i | i i i | [~ 

Standard 
(VA max > = 0.471 





I I I | I I l | I M | M l I | I I I | M 

AL 

<V/V mBX >=0.462 




_L 




± 



0.0 2.0 4.0 6.0 8.0 10.0 
-log(L/L Q ) 

Fig. 1. Luminosity function of halo white dwarfs for a stan- 
dard initial mass function (top panel) and a biased initial 
mass function (bottom panel). The observational luminos- 
ity function of halo white dwarfs is represented using a dot- 
ted line (Torres et al. 1998) and solid triangles, while the 
theoretical luminosity function is shown using a solid line 
and squares. See text for details. 



on the particular characteristics of the experiment and, 
hence, we consider different detection efficiencies and dif- 
ferent total exposures for the MACHO and EROS experi- 
ments. Specifically, for the case in which we analyze the re- 
sults of the MACHO collaboration we have taken 1.1 x 10 7 
stars during 5.7 yr and over 13.4 deg 2 , whereas the detec- 
tion efficiency has been modeled as: 



0.43e-( ln (*7T m )) 3 -/o.87 ; i>Tn 
£[ ) 1 o.43e-IM*7T m )| 234 /ii.i6 ; i<T„ 



(5) 



where T m = 250 days. This expression provides a good 
fit to the results of Alcock et al. (2000). For the EROS 
experiment we have used 0.7 x 10 7 stars over a wider field of 
84 deg 2 and over a period of 6.7 yr. Regarding the detection 
efficiency we have adopted a numerical fit to the results 
presented in Tisserand et al. (2006). 
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Fig. 2. Microlcnsing optical depth towards the LMC as a 
function of the limiting magnitude. Open and solid simbols 
represent the population of white dwarfs without and with 
the contribution of the ONe white dwarfs, respectively. The 
solid symbols have been shifted for the sake of clarity. 

3. Results 

3.1. The halo white dwarf luminosity function 

Despite the increasing number of surveys searching for 
white dwarfs - - like the Sloan Digital Sky Survey 
(Eisenstein et al. 2006), the 2 Micron All Sky Survey (Cutri 
et al. 2003), the SupcrCosmos Sky Survey (Hambly, Irwin & 
MacGillivray 2001), the 2dF QSO Redshift Survey (Vennes 
et al. 2002), and others — their success in finding halo 
white dwarfs has been limited. Thus, the observational de- 
termination of the halo white dwarf luminosity is still to- 
day rather uncertain. In fact, the two attempts to build 
such a luminosity function (Liebert et al. 1989; Torres et al. 
1998) have provided us only with the bright branch of the 
halo white dwarf luminosity function. Nevertheless, this is 
enough for our purposes, since we only need a normalization 
criterion and, hence, only an upper limit to the local den- 
sity of moderately bright dwarfs is needed. Consequently, 
we have used the luminosity function of Torres et al. (1998) 
and we have normalized the local density of white dwarfs 
obtained from our Monte Carlo simulations to its observed 
value, n ~ 9.0 x 10~ 6 pc~ 3 for log(L/i ) ^ -3.5 (Torres 
et al. 1998). 
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Table 1. Summary of the results obtained for the simulation of microlcnses towards the LMC for the MACHO model 
for an age of the halo of 14 Gyr, different model IMFs, and several magnitude cuts. 



From the distribution of white dwarfs obtained using 
our Monte Carlo simulations we compute the white dwarf 
luminosity using the 1/V max method (Schmidt 1968). It is 
important to mention that when deriving a luminosity func- 
tion using the 1/V max method a proper motion cut and a 
limiting magnitude are required. The set of selection cri- 
teria adopted here for computing the halo white dwarf lu- 
minosity function is the same as used in Garcfa-Berro et 
al. (2004). Namely, we have chosen a limiting magnitude 
my 1 = 17.5 mag and a proper motion cut /j, > 0.16" yr -1 . 
With all these inputs the luminosity functions in Fig. 1 
are obtained. The top panel shows the halo white dwarf 
luminosity function obtained using a standard initial mass 
function, whereas the bottom panel shows the luminosity 
function when the biased initial mass function of Adams 
& Laughlin (1996) is adopted. The simulated luminosity 
functions are represented as squares connected with solid 
lines, whereas the observational luminosity function is rep- 
resented as triangles connected with dashed lines. We also 
recall that, by construction, our samples are complete, al- 
though we only select about 10 white dwarfs using the se- 
lection criteria discussed before. However, our simulations 
do provide the whole population of white dwarfs, which is 
much larger. Hence, we can obtain the real luminosity func- 
tion by simply counting white dwarfs in the computational 
volume. This is done for all realizations and then we obtain 
the average. The result is depicted as a solid line in Fig. 
1. The true luminosity function steadily increases for lumi- 
nosities larger than \og{L/ Lq) ~ —5.0 and then sharply 
drops. This drop-off is given by the paucity of CO white 
dwarfs with appropriate ages (14 Gyr). Note however that 
the bulk of the population of ONe white dwarfs is located 
at much smaller luminosities, a consequence of the much 
shorter cooling timescales of these white dwarfs. In fact, 
for a typical halo age of 14 Gyr, the bulk of the ONe white 
dwarf population has already entered the fast Debye cooling 
phase and, consequently, would not be detectable with the 
current observational facilities. In the next sections we ex- 
plore if this elusive white dwarfs contribute significantly to 
the microlensing optical depth. It is also important to note 
that with the adopted limiting magnitude and proper mo- 
tion cut we obtain simulated white dwarf luminosity func- 
tions which are totally compatible with the observational 
one. Hence, the inclusion of massive ONe white dwarfs does 
not appreciably change the resulting white dwarf luminosity 
function, which is very similar to that obtained in Garcfa- 
Berro et al. (2004). 



3.2. Microlensing towards the LMC 

Firt of all, we analyze the result obtained by the MACHO 
collaboration. In Fig. 2 we show the contribution to the 
optical depth towards the LMC due to the white dwarf 
population as a function of the adopted limiting magni- 
tude. The results have been normalized to the value derived 
by Alcock et al. (2000), t = 1.2 x 10~ 7 . The open sym- 
bols represent the contribution if only CO white dwarfs 
are taken into account, while the solid symbols show the 
contribution to the microlensing optical depth when both 
CO and ONe white dwarfs are correctly included in the 
model white dwarf population. As can be seen, for none 
of the adopted initial mass functions the inclusion of the 
ONe white dwarf population does not significantly increase 
the contribution of white dwarfs to the microlensing optical 
depth towards the LMC, despite the fact that ONe white 
dwarfs arc much fainter than regular CO white dwarfs (see 
also Fig. 1). Specifically, the contribution of the white dwarf 
population is, respectively, of the order of 10% for the case 
of the standard initial mass function and somewhat larger 
(~ 15%) for the log normal initial mass function of Adams 
& Laughlin (1996). These figures are comparable to those 
already found in Garcfa-Berro et al. (2004). The only dif- 
ferences are that in the case of the standard mass function 
the contribution of ONe white dwarfs to the microlensing 
optical depth is clearly dominant only when the adopted 
limiting magnitude is of the order of 30, which is a totally 
unrealistic value. For the case of the log-normal initial mass 
function the results presented here show that the contribu- 
tion is nearly constant, independently of the adopted limit- 
ing magnitude, whereas when only the contribution of CO 
white dwarfs was considered the contribution to the opti- 
cal depth of the halo white dwarf population was clearly 
decreasing for increasing magnitude cuts. 

A summary of the results obtained with our Monte 
Carlo simulator can be found in Table 1, where we show 
for four selected magnitude cuts the number of microlens- 
ing events, the average mass of the microlenses, their av- 
erage proper motion, distance and tangential velocity, the 
corresponding Einstein crossing times and, finally, the con- 
tribution to the microlensing optical depth. It is important 
to discuss some of the numerical values in Table 1. For 
instance, it is clear that the larger the magnitude cut, the 
more massive the average mass of the lenses, as it should be 
expected from Fig. 1. In particular, for the case in which a 
standard initial mass function is used we obtain that for the 
largest limiting magnitude the average mass is ~ 0.9 Mq, 
indicating that in sizeable fraction of the Monte Carlo re- 
alizations the lens is an ONe white dwarf. Also, the log- 
normal initial mass function produces more microlensing 
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Fig. 3. Fraction of microlenses due to ONe white dwarfs 
with respect to the whole population of white dwarfs for 
the standard initial mass function — squares — and for 
the log-normal initial mass function of Adams & Laughlin 
(1996) — triangles. 



events, as one should expect, given that this biased initial 
mass funtion was tailored to produce more microlensing 
events. In fact for this initial mass function a maximum 
number of 6 microlensing events should be expected, while 
for the standard initial mass function we should expect 1 
microlensing event, at most. However, the contribution to 
the microlensing optical depth is only slightly larger for the 
Adams & Laughlin (1996) initial mass function. The reason 
for this is that the microlensing events for this distribution 
have shorter Eintein crossing times, as seen in table 1. 

The results obtained so far are not evident at first 
glance, since one may expect that ONe white dwarfs should 
be good microlensing candidates. As previously mentioned, 
ONe white dwarfs have a faster cooling rate than that of CO 
white dwarfs and, consequently, they reach much fainter 
magnitudes for the same cooling age. Hence, for reason- 
able halo ages one should naively expect that the probabil- 
ity that a ONe white dwarf could produce a microlensing 
event would be somewhat larger than that of a CO white 
dwarf, given that for reasonable halo ages practically all 
ONe white dwarfs have magnitudes larger than the mag- 
nitude cuts adopted here. However, even if this is indeed 
the case, we have shown that the total contribution of ONe 
white dwarfs is almost negligible. To clarify this result we 
have analyzed the fraction of microlenses due to ONe white 
dwarfs with respect to that of the total population. In Fig. 
3 wc show this fraction as a function of the limiting magni- 
tude for the two initial mass functions under study. As can 
be seen, the contribution of ONe white dwarfs is small for 
limiting magnitudes below 25 mag . Specifically, for the case 
of the standard initial mass function they only contribute 
a modest 2%, whereas for the log-normal initial mass func- 
tion the contribution is halved. This situation only reverses 
when magnitude cuts larger than ~ 27 mag are adopted. 
This result by itself is not explanatory of why the contribu- 
tion of ONe white dwarfs is not significant. We recall here 
that the contribution of an object to the total optical depth 
is given by Eq. (4), which depends on the Einstein crossing 
time which, in turn, depends on the Einstein radius and 
on the transverse velocity of the lens, £e = ^E/^tan- The 
Einstein radius scales as the root of the mass of the object 



and it also depends on the lens-object distance — see Eq. 
(3). We note that the average mass of an ONe white dwarf 
is larger than that of a CO white dwarf. Also, given the 
intrinsic faintness of ONe white dwarfs, their spatial dis- 
tribution in the computational volume is different because 
we are selecting microlensing candidates with magnitudes 
fainter than a given limiting magnitude. Thus, it can be 
expected that the contribution to the optical depth of an 
representative object of the these two populations should 
be different as well. 

In Table 2 we show the average parameters of the ONe 
white dwarf population susceptible to produce a microlens- 
ing event. The average mass of an ONe white dwarf is 
~ 1.1 M , while for a CO white dwarf it is ~ 0.6 M@. On 
the other hand, the average distance of ONe white dwarfs 
is in the range between about 2 and 4 kpc, independently of 
the limiting magnitude, while for the CO white dwarf pop- 
ulation the average distance increases for increasing magni- 
tude cuts. Finally, the average tangential velocities are very 
similar for all the magnitude cuts, given that the selection 
criteria are independent of the kinematical properties of the 
sample. With these data and using Eq. (3) and (4) the ratio 
of the contribution to the optical depth of a typical ONe 
white dwarf with respect to the contribution of a typical 
CO white dwarf is 



T~ONc _ toNc s(t 



CO) 



/ M ONc £>»£ c efco) 
tco too e(ioNo) V m co-Dql e(*ONc) 



(6) 



This ratio turns out to be roNe/Tco ~ 1-5- Recalling 
that the fraction r\ of ONe white dwarfs for limiting magni- 
tudes fainter than 25 mag is typically 0.02 for the standard 
initial mass function and 0.01 for the biased initial mass 
function, the increment in the total optical depth due to 
ONe white dwarfs can be estimated to be 



At 

to 



_ TQNc 
TCO 



(7) 



which represents an increment of roughly 3% for the case 
in which a standard initial mass function is considered and 
a 2% increment for the case of the log-normal initial mass 
function. These results are in nice agreement with those 
previously presented in Fig. 2. On the other hand, when 
the magnitude cut is 30 mag the fraction of ONc microlenses 
rj increases significantly and, thus, the fractional increase of 
the optical depth due to ONe white dwarfs consequently in- 
creases, reaching values as high as 100%. This fact is respon- 
sible for the different behaviour of the deepest magnitude 
bins of the left panel of Fig. 2, which show the situation 
for the standard initial mass function. The biased initial 
mass function suppresses the formation of moderately mas- 
sive ONe white dwarfs, and this is the reason why these 
faintest luminosity bins are not as populated as the equiv- 
alent bins for the case in which a standard mass function 
is considered. 

In a second set of Monte Carlo calculations we have 
simulated the observational data obtained by the EROS 
team. We recall here that the EROS collaboration have not 
found any microlensing event towards the LMC and one 
candidate event towards the SMC. Adopting a standard 
halo model and assuming tsmc = 1-4tlmc> the EROS re- 
sults imply an optical depth tq = 0.36 x 10 -7 (Tisserand et 
al. 2006), which is four times smaller than that obtained by 
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AL 






Magnitude 


17.5 


22.5 


27.5 


32.5 


17.5 


22.5 


27.5 


32.5 


(m) (M/M Q ) 


1.118 


1.106 


1.244 


1.130 


1.092 


1.082 


1.083 


1.101 


(d) (kpc) 


4.95 


3.98 


2.83 


3.80 


2.31 


2.02 


2.84 


5.75 


(Kan) (kmS- 1 ) 


253 


257 


250 


250 


266 


255 


250 


269 


(*e> (d) 


107.9 


91.6 


77.7 


104.1 


56.8 


61.6 


75.0 


99.0 



Table 2. Average values for the ONe white dwarf population. 
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AL 




Magnitude 


17.5 


22.5 


27.5 


17.5 


22.5 


27.5 


(Nwd) 


Oil 


Oil 


Oil 


1±2 


1±2 


0i2 


(m) (M/M ) 

M CV- 1 ) 


0.607 


0.595 


0.622 


0.631 


0.634 


0.642 


0.013 


0.011 


0.008 


0.034 


0.025 


0.010 


(d) (kpc) 


4.29 


4.52 


6.71 


1.50 


2.03 


5.39 


(Kan) (kmS" 1 ) 


256 


239 


246 


240 


244 


258 


(*b> (d) 


64.9 


77.0 


89.7 


37.9 


45.3 


74.8 


(r/r ) 


0.344 


0.372 


0.392 


0.368 


0.384 


0.505 



Table 3. Summary of the results obtained for the simulation of microlenses towards the LMC for the EROS model for 
an age of the halo of 14 Gyr, different model IMFs, and several magnitude cuts. 



the MACHO team. Although it is expected that the value 
of the optical depth obtained from our simulations should 
be only slightly different, it is as well true that this may 
be a test of the robustness of our numerical procedures. In 
particular, the detection efficiency of both experiments is 
very different. Additionally the areas (and the number of 
objects) surveyed by both teams are different. The data re- 
sults are summarized in table 3. Our simulations show that 
the white dwarf population could account for a 35% of the 
optical depth found by the EROS team if a standard ini- 
tial mass function is adopted, while for the non-standard 
initial mass function the contribution of the white dwarf 
population could be as large as 50%. On the other hand 
the expected number of objects has an upper limit of 1 
for the standard initial mass function and 2 for the log- 
normal initial mass function. Both results are in agreement 
with the results of the EROS experiment. Again, as it was 
the case for the simulation of the MACHO experiment, the 
contribution of ONe white dwarfs is small. All in all, it 
seems that the microlensing optical depth obtained by the 
MACHO collaboration is a clear overestimate. 

3.3. The Hubble Deep Field South 

Kilic et al. (2005) have recently re-observed the Hubble 
Deep field south (HDF-S), and have found three white 
dwarf candidates among several faint blue objects which 
exhibit significant proper motion and, thus, are assumed 
to belong to the thick-disk or the halo populations. 
Consequently, we have also performed a series of Monte 
Carlo simulations in the direction of the HDF-S (I = 
328.25° b = -49.21°) for a small window of 4.062 arcmin 2 . 
We have used the Johnson-Cousins U BVRI system instead 
of the WFPC2 photometry because the differences between 
both photometric systems is smaller than 0.02 mag for the 
range of colors under study (Holtzman et al. 1995). Also, no 
reddening was applied to the synthetic white dwarf stars. 
Contrary to what has been done until now the results pre- 
sented in this section are the average of 10 3 different re- 
alizations. Each of these realizations has been normalized 
to the local density of halo white dwarfs as previously de- 
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Fig. 4. Color-magnitude diagram for the white dwarf dis- 
tribution (ONe white dwarfs are circled) for the HDF-S 
of two typical simulations. The dashed line respresents the 
HDF-S observation limit. Also represented is the average 
expected location within lcr error of a typical ONe white 
dwarf. See text for details. 



8 



J. Camacho et al.: ONe white dwarfs and the Galactic halo 



scribed. The synthetic white dwarf population using this 
procedure is shown in the color-magnitude diagram of Fig. 
4. In this figure we represent two typical simulations for the 
halo white dwarf population in the direction of the HDF-S 
for the two initial mass functions under study. As can be 
seen, the number of white dwarfs susceptible to be detected 
in the HDF-S survey — that is, those with I magnitude 
smaller than 27 mag — is substantially larger for the log- 
normal initial mass function of Adams & Laughlin (1996) 
than for the standard initial mass function. Specifically, the 
average number of objects with / < 27 mag turns out to be 
6 ± 2 for the case in which a standar initial mass function 
is adopted, while for the log-normal initial mass function 
this number is 110 ±8. However, and in order to avoid con- 
fusion with blue extragalactic objects and main sequence 
stars, Kilic et al. (2005) restricted their search for white 
dwarfs candidates to colors in the range V—I < 0.4. Adding 
this new restriction we obtain that the expected number of 
white dwarfs should be 1 ± 1 for both initial mass functions. 
Although this result implies that both initial mass functions 
are compatible with the observations, the log-normal ini- 
tial mass function produces a large number of white dwarfs 
with colors in the interval 0.6 < V — I < 1.4, which has no 
observational counterpart. Additionally, in Fig. 4 we also 
show the only one ONe white dwarf obtained for each one 
of these two typical simulations. In both cases its location is 
shown as an encircled dot in the color-magnitude diagram. 
It is worth mentioning that in most of the 10 3 realizations 
an ONe white dwarf is found, and thus we also show the av- 
erage location of ONe white dwarfs in the color-magnitude 
diagram, along with the corresponding la error bars. Note 
that in any case ONe white dwarfs are much fainter and 
bluer than normal CO white dwarfs, as it should be ex- 
pected given that for a typical age of the halo most ONe 
white dwarfs have already reached the blue hook in the 
color-magnitude diagram. 

3.4. The dark matter density 

The results discussed so far indicate that, even in the case 
in which the contribution of ONe white dwarfs is taken 
into account, only a small fraction of the microlensing op- 
tical depth towards the LMC can be attributed to white 
dwarfs. We recall that if we adopt the microlensing opti- 
cal depth of the MACHO experiment this contribution is 
nearly a 20% for the biased initial mass function of Adams 
& Laughlin (1996) and <~ 10% for the standard initial mass 
function. Besides, for a spherical isothermal halo model the 
microlensing optical depth towards the LMC is given by 
the expression (Alcock et al 2000; Griest 1991): 

tlmc - 5.1 x 10- 7 / (8) 

where / is the fraction of the halo mass that is made of 
lensing objects. Thus, the white dwarf population would 
contribute / ~ 0.05 to the mass of the halo in the most 
optimistic case. 

However, we can go one step beyond using the results 
of our Monte Carlo simulations. In particular, we can com- 
pute the baryonic dark matter density in the form of white 
dwarfs using the l/V mSLX method. We proceed as follows. For 
each star of the sample we determine the maximum volume 





Standard 




AL 




CO ONe 


CO 


ONe 


ISM 


1.1 x W" 1 6.4 x 10"° 


2.8 x 10" 


' A 3.5 x 10~ D 


WD 


5.4 x 10" 5 9.5 x 10" 7 


9.2 x 10" 


4 5.3 x 10" 6 



Table 4. Density of baryonic matter (Mq/pc 3 ) in the 
Galactic halo within 300 pc from the Sun in the form gas 
returned to the interstellar medium (ISM) and in the form 
of white dwarfs (WD). 



over which each star can contribute as a microlensing event 
using the expression 

Lmax = "g"( r max — r min) (9) 

where r max is the radius of the volume in which we dis- 
tribute the objects of our sample, which in our case is the 
radius of Galactic halo, and r m ; n is the minimum volume 
for which a white dwarf still belongs to the sample consid- 
ering its apparent magnitude to be fainter than the adopted 
magnitude cut. Then, the number density of white dwarfs 
is 

» = E— ■ (10) 

i=1 ''max; 

Using this procedure we find that the contribution of 
white dwarfs to the baryonic dark matter would be roughly 
a 3% in the case in which a standard intial mass function is 
considered and nearly a 5% for the case in which the initial 
mass function of Adams & Laughlin (1996) is adopted. 

Finally, from our Monte Carlo simulations we can also 
derive the total density of baryonic matter in the Galactic 
halo within 300 pc from the Sun in the form of main se- 
quence stars, stellar remmants and in the corresponding 
ejected mass. We obtain p = 2.6 x 10~ 4 Mq P c ~ 3 for the 
standard initial mass function and 3.8 x 10 -3 Mq pc~ 3 for 
the log-normal intial mass function. The respective contri- 
butions of CO and ONe white dwarfs to the mass budget 
and of the mass returned to the interstellar medium are also 
shown in Table 4. Note that the total contribution of ONe 
white dwarfs is rather limited. The total density of baryonic 
matter obtained from our Monte Carlo simulations can be 
compared as well with the local dynamical matter density: 

v 2 t 

PDM = 47G%' (U) 

where v rot is the rotation velocity of the Galaxy and Rq is 
the Galactocentric distance. Thus, the fraction r\ of bary- 
onic matter of the Galaxy resulting from the white dwarf 
population can be estimated. Our results indicate that r\ 
would be a modest 0.02 for the case in which a standard 
initial mass function is adopted, whereas a sizeable frac- 
tion of the baryonic matter, rj = 0.52, can be accounted 
if the initial mass function of Adams & Laughlin (1996) is 
assumed. 

4. Conclusions 

We have analyzed the contribution of ONe white dwarfs 
to the MACHO content of the Galactic halo. We find that 
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although ONe white dwarfs fade to invisibility very rapidly 
and, thus, they are good baryonic dark matter candidates, 
their contribution to the microlensing optical depth towards 
the LMC is rather limited. In particular, we have found that 
when the contribution of ONe white dwarfs is taken into ac- 
count the microlensing optical depth does not increase sig- 
nificantly, independently of the adopted initial mass func- 
tion. If the microlensing optical depth is adopted to be that 
of the MACHO experiment, t = 1.2 x 10~ 7 (Alcock et al. 
2000) — which probably is an overestimate — we find that 
the fraction of the microlensing optical depth due to the 
whole white dwarf population is at most ~ 13% in the case 
in which a standard initial mass function is adopted and 
~ 19% if the log-normal initial mass function of Adams 
& Laughlin (1996) is considered. These values are roughly 
~ 3% larger than those already found by Garcfa-Berro 
et al. (2004), who only considered the contribution of CO 
white dwarfs. We have also studied if some of the candidate 
white dwarfs of the Hubble Deep Field South could be ONe 
white dwarfs and we have found that most probably this 
is not the case. Finally, we have also discussed the contri- 
bution of the whole white dwarf population to the mass of 
the Galactic halo. We have found that this contribution is 
of the order of a modest 5% in the most optimistic case. 
All in all, we conclude that white dwarfs are not significant 
contributors to the mass of the Galactic halo. 
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